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costing  the  U.S.  economy  over  S70  billion  annually.  Real-time  radiography  has 
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not  yet  demonstrated  a  significantly  better  ability  than  other  competing 
modalities  to  detect  the  early  stages  of  corrosion  when  corrective  actions  can 
most  efficiently  be  taken.  A  dual-energy  radiographic  technique  Is  reported 
which  will  significantly  enhance  the  corrosion  detection  and  damage  assessment 
capabilities  of  any  radiometric  system.  The  main  advantage  of  the  proposed 
method  Is  that  It  Improves  the  conspiculty  of  detail  In  a  digital  radiograph  by 
greatly  reducing  the  structural  noise  in  the  Image.  The  Improved  sensitivity 
offered  by  dual -energy  methods  can  be  traded  for  an  equivalent  detectability  at 
a  much  lower  dose  level,  an  Important  consideration  for  on-board  or  In-situ 
Inspections.  In  addition,  dual-energy  radiography  also  allows  the  user  to 
identify  materials  In  the  object  under  examination,  an  Important  consideration 
when  foreign  material  Is  present.  Results  presented  here  Indicate  that  the 
effective  atomic  number  can  be  determined  to  better  than  3X  absolute  error. 

This  newly  developed  dual-energy  technique  impacts  many  areas  of  corrosion 
inspection  and  diagnosis.  In  those  areas  where  real-time  radiography  has 
already  demonstrated,  or  Is  expected  to  demonstrate,  an  ability  to  produce 
useful  Inspections,  the  technique  will  serve  to  improve  the  sensitivity  of  the 
examination  and  perhaps  expand  the  range  of  valid  applications.  In  those 
Inspection  areas  where  real-time  radiography  offers  at  best  a  marginal  detection 
capability,  dual-energy  processing  may  make  some  otherwise  Impractical  applica¬ 
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detection  and  assessment  of  corrosion,  the  dual-energy  method  Is  considerably 
more  powerful  and  more  general  in  nature  and  could  rapidly  become  a  standard 
radiographic  procedure. 
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1.0  INTRODUCTION 


This  final  report  is  submitted  in  fulfillment  of  the  reporting  requirements 
of  AFWAL  Contract  No.  F33615-84-C-5089,  entitled  “A  New  Radiographic  Corrosion 
Inspection  Capability."  The  report  summarizes  the  results  of  an  ARACOR  SBIR 
Phase  I  effort  conducted  between  2  August  1984  and  8  April  1985  to  evaluate  the 
efficacy  and  sensitivity  of  a  dual-energy  radiographic  method.  A  two-pronged 
research  effort  was  developed  and  implemented,  providing  both  simulated  and 
actual  dual-energy  measurements.  Based  on  the  results  of  this  work,  it  is 
concluded  that  a  dual-energy  technique  can  significantly  enhance  the  corrosion 
detection  and  damage  assessment  capabilities  of  a  digital  radiographic  system. 

2.0  PROGRMf  DESCRIPTION 

Corrosion  of  equipment  and  material  is  a  ubiquitous  national  problem.  For 
the  Armed  Forces,  preventing,  detecting  and  repairing  corrosion  is  a  crucial 
part  of  their  continuous  struggle  to  keep  military  systems  in  full  combat 
readiness.  A  recent  study  by  the  National  Bureau  of  Standards,  "Economic 
Effects  of  Metallic  Corrosion  in  the  United  States,"  estimates  corrosion-related 
problems  cost  the  U.S.  economy  over  $70  billion  annually.  A  similar  figure  has 
not  been  placed  on  corrosion  losses  within  the  Armed  Forces,  but  it  is  certainly 
in  the  billions  of  dollars  (1).  In  October  1980,  the  Joint  Logistics  Commander 
established  a  Joint  Panel  on  Corrosion  Prevention  and  Control  with  one  of  its 
specific  mandates  being  to  provide  special  emphasis  on  the  development  of  new 
corrosion-sensitive  NOI  techniques  with  the  potential  to  provide  rapid, 
large-area  inspections.  Concurrently,  a  Joint  Technical  Coordinating  Group  on 
NDI  was  also  established  with  the  development  of  an  effective  corrosion 
inspection  technique  as  one  of  their  high  priority  objectives.  As  a  result  of 
these  and  other  related  efforts,  a  major  workshop  on  the  NDE  of  aircraft 
corrosion  was  held  at  Wright  Aeronautical  Laboratories,  Dayton,  Ohio,  in  May 
1983.  The  dual-energy  corrosion  detection  technique  reported  here  emerged  as  a 
direct  response  to  the  proceedings  of  that  workshop. 


2.1  Technical  Rationale 


Detection  of  corrosion  In  real  structures  is  complicated  by  the  myriad  of 
geometric  configurations  encountered  in  actual  objects  as  well  as  by  the  fact 
that  a  significant  amount  of  corrosion  is  either  not  immediately  accessible  to 
the  inspector  or  is  invisible  to  commonly  used  inspection  techniques. 
Consequently,  there  is  a  general  need  for  new  inspection  procedures  and/or 
technologies  which  are  more  quantitative,  more  reliable,  more  reproducible,  more 
rapid,  more  economical,  and  more  forgiving  than  current  methods.  Real-time 
radiography  has  been  identified  as  one  of  several  technologies  that  offers 
significant  promise  of  providing  improved  corrosion  detection  (1-4).  Real-time 
systems  have  the  potential  to  provide  rapid,  semi-quantitative  inspection;  and 
if  the  data  are  converted  to  digital  format  at  some  stage  in  the  signal 
processing  chain,  various  spatial-  and  temporal -mask  subtraction  techniques  and 
a  variety  of  enhancement  routines  could  conceivably  be  used  to  provide  an 
automatic  corrosion  detection  and  assessment  capability.  Yet,  in  spite  of  the 
potential  advantages  allegedly  possessed  by  real  time  systems,  they  have  not  to 
date  demonstrated  a  significantly  better  ability  to  detect  corrosion  in  its 
early  stages  (when  corrective  actions  can  most  efficiently  be  taken)  than  other 
competing  modalities  (2).  ARACOR,  which  had  been  exploring  dual-energy  imaging 
methods  as  they  apply  to  Computed  Tomography  for  more  than  half  a  decade, 
recognized  the  possibility  of  significantly  improving  corrosion  detection  and 
assessment  technology  by  applying  dual-energy  techniques  to  digital  radiographic 
systems. 

2.2  PrograM  Objectives 

The  overall  objective  of  the  work  reported  here  was  to  provide  both 
simulated  and  actual  demonstrations  of  dual-energy  radiographic  corrosion 
inspections.  Software  simulations  are  the  most  efficient  means  of  rapidly 
ascertaining  the  general  sensitivity  of  a  new  inspection  method,  and  actual 
measurements  of  real  components  are  the  most  definitive  means  of  unequivocally 
establishing  feasibility.  As  a  result,  a  two-pronged  research  approach  was 
conceptualized  and  implemented.  The  specific  task-oriented  objectives  which 
were  established  and  the  rationale  behind  them  can  be  enumerated  as  follows: 
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*  Perfora  Coaputer  Siaulations  of  Dual-Energy  Corrosion  Detection 


.  Most  efficient  method  of  determining  the  general  sensitivity  of 
the  dual-energy  technique. 

.  Crucial  for  guiding  actual  experiments  in  the  most  fruitful  and 
useful  directions. 

.  Required  to  provide  an  understanding  of  the  computer  architecture 
and  software  requirements  that  will  be  needed  in  an  actual  system. 

*  Perform  Actual  Dual -Energy  Measurements  on  a  Suitable  System 

.  Vital  for  corroborating  the  results  of  the  computer  simulations 
and  for  validating  the  parametric  studies  performed. 

.  Necessary  to  provide  unequivocal  proof-of-principle  demonstration 
of  dual-energy  methods. 

.  Needed  to  uncover  unanticipated  problem  areas  so  that  a  practical 
design  of  a  dual-energy  system  can  be  more  efficiently  evolved. 

The  program  effort  was  successful  in  accomplishing  all  of  the  goals  set  for 
it.  The  results  of  the  software  simulations  are  presented  in  the  next  section 
and  the  results  of  the  actual  measurements  in  Section  4.  A  summary  of  the 
conclusions  and  recommendations  reached  by  this  effort  are  given  in  Section  5. 

2.3  Technical  Overview 

The  term  "dual -energy"  is  utilized  in  this  context  to  indicate  that  as  a 
standard  part  of  the  inspection  procedure,  radiographs  are  taken  at  two 
different  effective  energies  and  implicitly  assumes  that  the  radiographs  are 
digital  at  some  point  in  the  processing  so  that  they  may  be  analyzed  to  provide 
additional  information  not  otherwise  available  with  conventional  (single-energy) 
radiographic  systems.  The  approach  relies  on  the  fact  that  in  the  energy  range 
from  a  few  keV  to  many  MeV,  all  but  a  small  number  of  minor  effects  are 
explainable  in  terms  of  just  three  photon-matter  interactions  --  photoelectric, 
Compton,  and  pair  production;  and  that  at  energies  below  about  1  MeV,  the 
Compton  effect  and  the  photoelectric  effect  dominate,  while  above  several  MeV, 
the  Compton  effect  and  pair  production  dominate'*’.  Consequently,  since  no  more 


"^An  excellent  introductory  explanation  of  the  physics  of  these  processes  can 
be  found  in  (5)  and  a  more  advanced  treatment  in  (6). 


than  two  interactions  need  to  be  considered  at  one  time  to  account  for  the 
attenuation  of  x-rays  in  matter  when  dealing  with  the  physics  of  radiography, 
digital  radiographs  obtained  at  two  different  energies  can  be  analyzed  to  obtain 
the  individual  contributions  of  the  competing  interactions  to  the  total  attenu¬ 
ation  (7).  Further,  because  the  photoelectric  and  pair  production  effects  are 
particularly  sensitive  to  atomic  number  and  the  Compton  effect  is  particularly 
sensitive  to  density,  the  effective  atomic  number  and  density  of  each  point  in  a 
radiograph  can  be  determined  from  projections  taken  at  two  different  energies. 

In  other  words,  detailed  material  properties  can  be  extracted  from  dual -energy 
radiographs. 

To  get  a  sense  of  how  this  is  accomplished,  consider  in  somewhat  more 
detail  the  physical  basis  for  dual-energy  digital  radiography.  Broadly  speak¬ 
ing,  the  dual-energy  technique  relies  upon  two  key  approximations  to  the  physics 
of  the  x-ray  attenuation  process.  First,  the  assumption  is  made  that  the  local 
mass  absorption  coefficient  of  the  absorbing  material,  /x,  can  be  adequately 
described  in  terms  of  only  two  attenuation  mechanisms,  and  i.e., 

/x{E,Z.A)  «Mi(E,Z,A)+  /^2(E.Z»A)  , 

where  the  full  dependence  of  n  on  the  x-ray  energy,  E,  and  the  local  atomic 
number,  Z,  and  the  local  mass  number.  A,  of  the  material  has  been  explicitly 
displayed.  As  stated  above,  describing  as  the  sum  of  only  two  attenuation 
mechanisms  is  a  good  approximation  for  most  radiographic  situations.  Second, 
the  assumption  is  made  that  each  attenuation  coefficient  can  be  adequately 
described  as  the  product  of  two  terms,  one  of  which  depends  solely  on  the  energy 
of  the  incident  photon,  the  other  of  which  depends  strictly  on  the  atomic  number 
of  the  interacting  matter,  and  neither  of  which  carries  an  A  dependence,  i.e., 

/Xi(E,Z,A)  Mi(E,Z)  T?i(E)  ai(Z)  , 
where  V  represents  the  energy-dependent  term  and  a  the  Z-dependent  term. 


The  dual-energy  method  leverages  these  approximations  by  recognizing  that 
the  description  of  the  x-ray  transport  of  a  polychromatic  beam  through  matter 
involves  both  an  average  over  the  energies  in  the  incident  flux  and  an  integral 
along  the  ray  path.  Since  the  E  and  Z  dependencies  separate  by  construction, 
the  line  integral  of  along  the  beam  path  can  be  written  as: 

/At(E,Z)ds  =T?i(§)/ai(Z)ds  +»?2(^)/®2(Z)ds  . 

where  ds  is  an  element  of  length  along  the  ray  path  of  the  beam  and  ^i(£)  and 
^2(^)  are  two  constants  that  depend  only  on  the  effective  energy,  E,  of  the 
fluence.  The  value  of  this  fact  is  that  7?x  and  172  can  be  determined  empirically 
by  suitable  calibration  procedures  and  used  to  algebraically  process  two 
radiographs  taken  at  different  x-ray  tube  energies  to  extract  a  pair  of  new 
images,  one  that  depends  on  folds  and  the  other  on  J‘a2ds,  Since  01  and  a2  depend 
only  on  the  material  properties  of  the  object  of  interest,  they  provide  a  method 
of  determining  the  additional  information  about  the  material  under  examination. 

One  particularly  powerful  use  of  this  technique  is  to  “cancel"  the 
predominant  material  in  the  image  to  improve  the  conspicuity  of  other  materials 
that  may  be  present.  This  is  accomplished  by  performing  a  weighted  subtraction 
of  the  01  image  from  the  02  image,  the  relative  proportions  being  selected  so 
that  oi  =  02  for  the  material  the  user  wishes  to  eliminate.  Since  the  magnitude 
of  01  and  02  vary  systematically  over  the  elements,  in  general  only  one  element 
(or  compound)  at  a  time  will  be  cancelled  with  this  technique.  Materials  in  the 
subtracted  image  with  a  higher  atomic  number  will  have  a  negative  residual; 
those  with  a  lower  atomic  number  will  have  a  positive  residual  (or  vice  versa, 
since  the  sign  is  arbitrary). 

2.4  Significance  of  the  Dual-Energy  Technique 

The  power  of  the  dual-energy  technique  is  twofold:  (1)  it  improves 
conspicuity  by  eliminating  structural  noise  in  the  image;  and  (2)  it  allows  the 
user  to  identify  different  materials  in  the  object  under  examination.  For 
instance,  if  an  aircraft  structure  is  composed  largely  of  aluminum,  then  the 


aluminum  can  be  cancelled  and  the  residual  image  studied.  Only  non-aluminum 
material  would  be  present  in  the  new  image.  Normally  expected  items,  such  as 
fasteners,  gaskets,  sealants,  etc.,  would  be  readily  recognized  by  the 
inspector;  and  any  corrosion  or  other  foreign  matter,  if  it  exists,  would  be 
visible  as  anomalous  patches  of  unexpected  material.  Since  non-aluminum  items 
are  generally  well  known  and  well  specified;  and  since  further,  those  areas  most 
likely  to  corrode  are  presumably  well  defined,  it  should  be  a  straightforward 
task  to  search  for  the  presence  of  corrosion  in  susceptible  components.  This 
technique  might  also  turn  out  to  offer  a  rapid  way  to  screen  for  missing 
components,  such  as  critical  washers  or  gaskets,  since  the  conspicuity  of  small 
parts  is  greatly  increased  when  more  dominant  materials  are  cancelled  in  a 
radiograph. 

Further,  dual-energy  radiography  can  also  be  used  in  a  more  general  way  to 
identify  materials.  Starting  with  a  normal  image,  the  user  can  sequentially 
cancel  one  material  after  another  while  watching  the  residual  image.  When  a 
particular  structure  or  component  suddenly  disappears  from  the  image,  the 
operator  can  infer  its  material  composition  from  the  particular  weighting  used 
to  achieve  the  cancellation.  Although  this  may  not  be  especially  useful  for 
identifying  standard  components  where  a'  priori  information  can  more  readily  be 
used  to  Identify  material  composition,  it  could  on  the  other  hand  be  very 
important  in  identifying  foreign  substances,  differentiating  contaminants,  such 
as  water  from  fuel,  or  in  assessing  the  degree  of  corrosion  present.  The 
results  of  this  work  indicate  that  the  effective  atomic  number  of  materials  can 
be  determined  to  better  than  three  percent  absolute  error. 

Thus,  the  dual-energy  technique  offers  a  powerful  means  of  dealing  with  a 
large  subset  of  corrosion  detection  problems.  Foremost,  it  should  significantly 
enhance  or  augment  the  corrosion  detection  and  damage  assessment  capabilities  of 
any  real-time  system.  Furthermore,  the  improved  sensitivity  offered  by 
dual-energy  methods  can  be  traded  for  an  equivalent  detectability  at  a  much 
reduced  dose  level,  an  important  consideration  in  close  environments,  such  as 
on-board  or  in-situ  inspections.  In  those  areas  where  real-time  radiography  has 
already  demonstrated,  or  is  expected  to  demonstrate,  an  ability  to  produce 
useful  inspections,  the  proposed  technique  will  serve  to  improve  the  examination 
and  to  expand  the  range  of  valid  applications.  In  those  areas  where  real-time 
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radiography  is  at  best  marginal,  dual-energy  processing  may  make  some  otherwise 
impractical  applications  feasible.  Either  way,  wherever  real-time  radiography 
finds  an  application  in  corrosion  detection  and  prevention,  dual-energy 
radiography  will  offer  substantial  improvements  in  detection  sensitivity  and 
damage  assessment. 

3.0  SOFTIMRE  DEVELOPNENT 

The  goal  of  the  software  development  effort  was  to  generate  a  set  of 
algorithms  which  would  allow  digital  radiographs  of  objects  taken  with  two 
distinctly  different  bremsstrahlung  spectra  to  be  decomposed  into  images  based 
on  the  competing  attenuation  mechanisms.  A  secondary  goal  was  to  perform 
detailed  simulations  to  aid  the  development  and  provide  an  initial  indication  of 
the  sensitivity  of  the  technique. 

The  major  tasks  needed  to  achieve  these  goals  were: 

*  Selection  of  the  parameterization  to  be  used  in  the  image 
decomposition; 

*  Generation  of  an  analytical  representation  of  the  calibration 
radiographs  in  terms  of  the  selected  parameters;  and 

*  Development  of  appropriate  linear  combinations  of  the  parameter 
images  to  allow  elemental  cancellations. 

The  results  of  these  subtasks  are  presented  in  the  following  subsections. 

3.1  Selection  of  Parameters  for  Image  Decomposition 

The  concept  of  elemental  analysis  by  dual-energy  radiography  originated  in 
the  medical  radiology  community(8) .  The  x-ray  attenuation  process  in  biological 
specimens  at  medical  x-ray  energies  is  dominated  by  Compton  and  photoelectric 
processes.  Radiologists  realized  that  since  the  Compton  cross-section  is  nearly 
independent  of  x-ray  energy  (E)  and  atomic  number  (Z)  while  the  photoelectric 
cross-section  varies  approximately  as  E-3  and  Z"3,  it  should  be  possible  to 
distinguish  the  two  attenuation  processes  from  a  pair  of  radiographs  taken  at 
different  energies.  However,  early  attempts  did  not  produce  particularly  good 
results  for  two  reasons:  (1)  initial  models  of  the  attenuation  mechanism  were 
not  adequate  to  handle  the  presence  of  a  weaker,  but  non-negligible,  third 


attenuation  mechanism  (I.e.,  Rayleigh,  or  coherent,  scattering);  and  (2)  the 
numerical  methods  of  processing  dual-energy  data  had  not  yet  evolved  to  the 
point  where  results  could  be  obtained  efficiently.  With  the  subsequent 
development  of  more  efficient  decomposition  schemes  and  data  processing 
techniques,  medical  dual-energy  methods  eventually  evolved  workable  approaches. 
While  these  techniques  produced  adequate  results  with  biological  materials  and 
relatively  low-energy  x-ray  sources  characteristic  of  medical  systems,  they  were 
Inadequate  for  Industrial  or  military  applications  In  which  the  predominant 
materials  consist  of  hIgher-Z  elements  and  radiometric  densities  require 
higher-energy  sources  to  obtain  adequate  signal-to-nolse  levels. 

For  the  range  of  materials  of  Interest  here,  the  absorption  process  Is  fed 
by  three  Interactions,  the  Compton,  the  coherent  (Rayleigh)  and  the  photoelec¬ 
tric.  The  Z  and  E  dependences  of  all  three  of  these  components  are  complex. 

The  coherent  cross-section  is  comparable  to  the  Compton  at  low-energy  and  com¬ 
parable  to  the  photoelectric  at  high-energies.  Just  as  It  did  for  early  medical 
attempts,  the  existence  of  a  third  Interaction  complicates  matters  significantly 
If  not  treated  properly.  The  approach  taken  here  was  to  model  the  attenuation 
process  In  terms  of  two  fictitious  cross  sections  which  give  a  good  fit  to  the 
mass  absorption  coefficients  without  trying  to  make  a  firm  connection  to  the 
underlying  details  of  the  processes  Involved.  The  terminology  Compton  and 
photoelectric  has  been  retained,  however,  because  at  high  energies  the  one  cross 
section  Is  essentially  the  same  as  the  Compton  Interaction  while  at  low  energies 
the  other  cross  section  Is  essentially  the  same  as  the  photoelectric 
interaction. 

3.2  Development  of  Calibration  Procedures 

Once  a  parameterization  is  selected.  It  Is  necessary  to  calibrate  the 
parameterization  by  scanning  an  object  of  known  composition  and  thickness.  The 
approach  used  here  was  to  describe  the  response  of  the  system  to  the  known  cali¬ 
bration  object  In  terms  of  the  established  model  and  thereby  deduce  the  various 
free  parameters  in  the  model.  Step  wedges  of  titanium  and  graphite  were  selec¬ 
ted  for  the  calibration  phantom  and  oriented  perpendicular  to  each  other.  Grap¬ 
hite  was  chosen  because  the  photoelectric  contribution  to  Its  x-ray  absorption 
Is  negligible  above  40  keV.  The  choice  of  titanium  was  a  compromise  between  the 
practical  requirements  for  a  material  whose  absorption  is  predominantly 
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photoelectric  and  whose  mechanical  properties  permit  It  to  be  readily  machined. 
Additionally,  the  two  materials  bracket  the  range  of  atomic  numbers  of  Interest 
and  are  also  commercially  available  In  highly  homogeneous  forms  with  excellent 
batch-to-batch  uniformity. 

Simulated  high-  and  low-energy  Images  of  the  calibration  phantom  were  then 
generated.  The  simulation  used  a  calculated  bremsstrahlung  spectrum.  Including 
the  tungsten  characteristic  line  radiation  (k^  and  k^),  based  on  the  model  of 
Reference  9.  The  spectrum  was  calculated  at  1  keV  Intervals  from  20  keV  to  the 
assumed  tube  voltage  (250  keV  for  the  low-energy  spectrum,  and  420  keV  for  the 
high-energy  spectrum).  The  high  energy  spectrum  was  presumed  to  be  filtered 
with  3  mm  of  brass  to  Improve  the  contrast  between  the  two  Images.  The  x-ray 
transmission  was  calculated  for  each  energy  bln  and  the  results  summed  to  give 
an  Integrated  transmitted  Intensity  for  each  graphite-titanium  combination.  The 
transmission  calculations  were  performed  using  the  mass  absorption  coefficients 
from  Reference  10.  Photon  noise  was  also  added  to  each  measurement  assuming 
Poisson  statistics. 

The  assumption  was  then  made  that  the  Compton  and  photoelectric  components 
of  the  low-energy  and  high-energy  Images  of  the  calibration  phantom,  C  and  P 
respectively,  could  be  represented  to  sufficient  accuracy  In  the  form: 

C(x,y)  »  fl(H(x,y),L(x,y)) 

P(x,y)  =  f2(H(x,y),L(x,y)), 


.v’.i 


where  fl  and  f2  are  polynomials  In  H(x,y)  and  L(x,y),  the  normalized  high-  and 
low-energy  Images,  respectively.  The  advantage  of  this  formulation  Is  that  If 
suitable  functions,  fl  and  f2,  can  be  found,  an  especially  rapid  method  of 
decomposing  measurements  Into  their  component  Images  can  be  employed.  Two 
orthogonal  functions,  fl  and  f2,  were  chosen  for  their  computational  properties 
and  a  code  was  developed  to  extract  the  coefficients  of  the  functions  from  the 
calibration  data.  The  simulated  calibration  Images  were  submitted  to  this 
calibration  routine  and  the  results  were  then  evaluated  by  generating  Compton 
and  photoelectric  Images  for  a  number  of  test  cases.  Some  results  are  presented 
in  Section  3.4. 
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3.3  Construction  of  Residual  loages 

The  Compton  and  photoelectric  images  extracted  from  the  simulated 
radiographs  were  further  processed  to  demonstrate  the  cancellation  of  a 
particular  material  In  an  image.  This  was  accomplished  by  talcing,  on  a 
point-by-point  basis,  a  linear  combination  of  the  component  images.  Thus,  if 
C(Z)  is  the  Compton  attenuation  of  a  particular  material  and  C(Z)  =  l((Z)  P(Z), 
where  P(Z)  is  the  photoelectric  attenuation  and  k(Z)  is  a  constant  whose  value 
is  specific  to  the  material  to  be  cancelled,  then  the  residual  at  each  point  in 
the  image  can  be  calculated  as  follows: 

R  «  C(Z)  -  lc(Z)P(Z). 

With  this  prescription,  the  material  of  interest  will  vanish.  If  the  object 
contains  material  with  a  lower  atomic  number  than  Z,  a  positive  residual  will 
result.  If  it  contains  material  with  a  higher  atomic  number,  the  residual  will 
be  negative.  This  is  illustrated  in  the  examples  given  in  Section  3.4. 

A  related  approach  is  to  form  the  best  visual  cancellation  and  then  obtain 
an  estimate  of  Z  using  parametric  forms  for  C(Z)  and  P(Z).  An  analysis  of  the 
inherent  errors  in  this  approach  reveal  that  the  uncertainties  are  dominated  by 
the  uncertainty  in  the  assumed  value  of  A  for  a  given  Z.  A  review  of  the  low-Z 
elements  showed  that  the  average  value  of  Z/A  for  the  naturally  occurring  stable 
elements  is  0.474  tO.021,  for  2«Z«30.  There  are  several  competing  mechanisms 
which  cause  <Z/A>  to  fluctuate  to  this  degree:  (1)  hydrogen  has  an  anomalously 
high  value  of  Z/A  »  1,  which  drives  <Z/A>  for  compounds  with  large  hydrogen 
content  strongly  upwards;  (2)  low-Z  compounds  (i.e.,  Z'^20,  which  includes  the 
corrosion-related  compounds  as  well  as  the  organics  and  many  inorganics)  are 
dominated  by  elements  with  Z/A  ~0.5;  and  (3)  higher-Z  compounds  (i.e., 
20<Z<30,  which  include  the  commonly  occurring  iron-series  elements:  iron, 
copper,  and  zinc)  are  dominated  by  elements  for  which  Z/A  ~0.46.  Other  sources 
of  error  include  the  accuracy  of  the  parameterization  (~1X)  and  the  ability  of 
the  observer  to  empirically  adjust  the  ratio  of  Compton  and  photoelectric  images 
(  ~Ul).  When  all  sources  of  error  are  analyzed  and  quantitatively  folded 
together,  the  ability  of  dual-energy  radiography  to  identify  an  unknown  material 
is  found  to  be  roughly  three  percent. 


3.4  Presentation  of  Slaulatlon  Results 


Simulated  high-energy  and  low-energy  images  of  a  16  step  x  16  step 
calibration  phantom  were  generated.  These  images  are  shown  in  Figure  1.  The 
graphite  thickness  is  constant  horizontally  with  the  thinnest  step  at  the  bottom 
of  the  image.  The  titanium  thickness  is  constant  vertically  with  the  thinnest 
step  at  the  left. 

The  coefficients  of  the  calibration  functions  were  extracted  and  used  to 
generate  Compton  and  photoelectric  images  of  the  step  wedge  phantom.^  These 
images  are  shown  in  Figure  2.  Because  the  photoelectric  cross  section  of 
graphite  is  negligible,  the  photoelectric  image  is  a  good  representation  of  the 
titanium  component  of  the  phantom.  On  the  other  hand,  since  the  titanium  has  a 
significant  Compton  component  in  its  absorption,  the  Compton  image  looks  very 
similar  to  the  high-energy  image.  However,  if  an  appropriate  fraction  of  the 
photoelectric  image  is  subtracted  point-by-point  from  the  Compton  image,  the 
effect  of  the  titanium  can  be  removed  from  that  image. 

In  the  units  used  for  the  simulation,  C(Ti)  *  .6988  *  P(Ti).  By  forming 
the  following  linear  combination  of  images,  a  residual  image  which  features  the 
graphite  is  obtained,  i.e.: 

R(graphite)  *  C(x,y)  -  .6988  *  P(x,y), 

where  R  stands  for  the  "residual"  image.  The  result  of  the  subtraction  is  shown 
in  the  left-hand  portion  of  Figure  3.  A  similar  procedure  can  be  used  to 
eliminate  the  graphite  and  feature  the  titanium.  This  is  shown  in  the 
right-hand  side  of  Figure  3.  The  quality  of  the  subtractions  is  an  indication 
of  the  validity  of  the  approach. 


^This  is  not  a  case  of  "data  in,  data  out."  The  simulated  images  were 
generated  from  information  in  References  9  and  10  and,  therefore,  totally 
independent  of  the  parameterization  selected  here  on  the  basis  of  convenience, 
functionality,  and  performance. 
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Figure  2.  Compton-photoelectric  separation  of  simulated  calibration 
phantom  radiographs. 
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To  test  the  sensitivitjrof  the  duiT-energy  Imaging  technique,  another  test 
was  performed,  this  time  using  two  new  crossed  step  wedges,  one  aluminum  and 
silicon  and  the  other  aluminum  and  aluminum-hydroxide.  The  Images  were  combined 
In  a  single  simulation  by  discarding  every  other  aluminum  step,  producing  a 
phantom  whose  lower  half  represents  crossed  steps  of  aluminum  and  aluminum 
hydroxide  and  the  upper  half  crossed  steps  of  aluminum  and  silicon.  Figure  4 
shows  the  simulated  high-  and  low-energy  radiographic  Images  with  noise  added. 
Figure  5  shows  the  Compton  and  photoelectric  Images  extracted  from  the  high-  and 
low-energy  Images.  Because  all  three  materials  have  a  significant  photoelectric 
component,  these  Images  differ  little  from  the  original  radiographs.  From  a  fit 
of  the  aluminum  mass  absorption  coefficient.  It  can  be  shown  that  C(A1)  ^  4.376 
*  P(A1).  Thus  the  residual  Image  formed  by  the  poInt-by-poInt  subtraction, 

R  -  C(x.y)  -  4.376  *  P(x,y), 

would  be  expected  to  cancel  the  aluminum  In  the  Image  and  obtain  a  positive 
residual  If  the  object  contains  material  with  a  lower  Z  than  aluminum  and  a 
negative  residual  If  the  object  contains  material  with  a  higher  Z  than 
aluminum.  The  resultant  positive  and  negative  residual  Images  are  shown  In 
Figure  6.  With  the  aluminum  cancelled,  the  silicon  and  aluminum  hydroxide  step 
wedges  are  clearly  resolved. 

Next,  .aaJoULIaL measure -of  the  sensitivity  of  the  technique  was  obtained 
by  empirically  adjusting  the  relative  proportions  of  the  Compton  and  photo¬ 
electric  Images  required  to  cancel  first  the  silicon  step  wedge  and  then  the 
aluminum  hydroxide  one  (see  Figure  7).  The  ratios  corresponding  to  the  best 
visual  cancellation  were  then  Inverted  to  obtain  an  estimate  of  Z. 

The  empirically  adjusted  ratio  for  a  silicon  cancellation  was  3.503  and  for 
aluminum  hydroxide  7.548.  These  ratios  produced  a  calculated  value  of  Z  * 
14.210.2  for  silicon  an'd  ifV4l0.2  for  aluminum  hydroxide.  Both  agree  within 
the  expected  error  with  the  accepted  value  of  Z  (or  In  the  latter  case,  the 
effective  value  of  Z). 


ICO 

gra 


A 


y. 

>, 

I 


I 


I 


This  appears  to  be  a  good  estimate  of  the  typical  sensitivity  of  the 
technique  in  a  practical  situation  when  no  other  information  is  utilized.  On 
the  other  hand,  if  additional  a  priori  information  is  available,  further 
improvement  is  possible.  One  approach  is  to  use  an  iterative  calculation  to 
obtain  a  better  estimate  of  Z.  To  do  this,  a  higher  value  of  <Z/A>  is  used  when 
the  initial  estimate  of  Z  is  less  than  20'*'  and  a  lower  value  if  Z  is  greater 
than  20.  As  an  example,  when  the  value  of  Z  ^  14.2  was  obtained,  a  more 
reasonable  value  of  <Z/A>  =  0.5  was  used  to  recompute  a  value  of  Z  on  the 
grounds  that  most  elements  with  a  Z  on  the  order  of  14  tend  to  have  a  value  of 
Z/A  nearer  0.5.  The  once-iterated  value  of  Z  obtained  in  this  way  was  14.0,  the 
correct  value  for  silicon.  The  same  calculation  for  the  aluminum  hydroxide  gave 
a  revised  estimate  of  Z-effective  of  11.2,  which  compares  acceptably  with  the 
"true"  value  of  10.8. 
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^Roughly  speaking,  <Z/A>  equals  the  mean  value  of  0.474  given  above  for  values 
of  Z  near  20.  (See  Section  3.3). 
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4.1  Overvie* 


To  test  the  procedures  developed  with  the  aid  of  simulated  data  on  real 
images,  access  to  a  digital  radiographic  system  was  arranged.  The  Air  Force 
Advanced  Computed  Tomography  System  Number  One  (AF/ACTS-I)  was  selected  on  the 
basis  of  its  proximity  and  availability  to  ARACOR,  our  familiarity  with  its  data 
acquisition  system,  and  the  unit's  ability  to  produce  high  quality  dual-energy 
radiographic  images.  The  dual-energy  methods  developed  here,  however,  are  quite 
general  and  the  demonstration  could  have  been  conducted  on  any  real-time,  or 
near-real-time,  system  capable  of  producing  digital  radiographs  at  different 
effective  energies. 

Several  typical  specimens  of  corroded  parts  were  provided  for  study  by 
A1  Rogel,  Nondestructive  Inspection  Manager,  McClellan  AFB.  The  item  selected 
for  the  most  detailed  testing  was  an  aluminum  antenna  panel  with  subsurface 
corrosion  within  stress- induced  delaminations.  To  complement  this  object, 

ARACOR  also  developed  a  special  test  phantom  for  analyzing  sensitivity.  Results 
of  these  studies  are  presented  in  this  section. 

4.2  Calibration  Measurements 

The  calibration  phantom  consisted  of  crossed  step  wedges,  one  of  graphite 
and  the  other  of  titanium,  each  comprised  of  16  steps.  The  titanium  steps  were 
10.2-cm  long  by  0.625-cm  wide  while  the  graphite  steps  were  10.80-cm  long  by 
0.597-cm  wide.  Thus  there  was  a  segment  of  each  graphite  step  which  was  not 
overlaid  with  titanium,  and  a  segment  of  each  titanium  step  which  was  not 
overlaid  with  graphite.  This  construction  resulted  in  a  total  of  17  x  17 
distinct  titanium-graphite  thickness  combinations  or  elements  (see  Figure  8). 

The  calibration  measurements  were  conducted  by  highly  oversampling  each 
element  of  the  calibration  phantom  at  both  energies.  The  low-energy  scans  were 
performed  at  250  kV  with  no  additional  filtration;  the  high-energy  ones  at 
420  kV  with  3-mm  brass  additional  filtration.  The  data  for  each  "element"  of 


the  phantom  were  then  averaged  before  being  Input  to  the  calibration  routine. 
Because  of  finite  beam  size  and  the  fact  that  the  beam  traversed  the  step  wedges 
at  a  different  angle  for  each  detector  array  element,  it  was  necessary  to 
develop  an  algorithm  to  select  only  the  data  points  for  which  the  entire  beam 
passed  through  the  full  thickness  of  each  graphite-titanium  element  and  to  apply 
a  correction  for  the  apparent  thickness  due  to  the  different  angles  of 
incidence.  The  data  set  passed  to  the  calibration  program  consisted  of  an 
average  of  a  minimum  of  500  data  points  for  the  thickest  steps  and  up  to  a 
maximum  of  1100  data  points  for  the  thinnest  steps.  The  standard  deviations  of 
the  mean  values  of  the  data  were  typically  less  than  one  percent,  and  the  fits 
obtained  from  the  calibration  code  were  found  to  be  excellent. 

4.3  Test  Object  Scans 

Radiographic  image  data  were  obtained  for  a  number  of  objects.  Three 
examples  are  shown  here. 

4.3.1  Dual-Energy  Calibration  Phantom.  Figure  9  shows  high-  and  low- 
energy  digital  radiographs  of  the  calibration  phantom  taken  with  the  tube  oper¬ 
ated  at  250  and  420  kV.  Figure  10  shows  the  extracted  Compton  and  photoelectric 
images.  As  with  the  simulated  images,  the  titanium  step  wedge  is  well  defined 
in  the  photoelectric  image,  whereas  the  Compton  image  is  not  materially  differ¬ 
ent  from  the  original  images.  Figure  11  shows  the  result  of  subtracting  0.7 
times  the  photoelectric  image  from  the  Compton  image.  As  in  the  simulations, 
the  graphite  step  wedge  is  now  well  developed  in  the  resulting  image. 

4.3.2  Alimilnua— Sapphire — Silicon  Composite  Phantom.  A  pair  of  aluminum 
step  wedges  used  tn  some  preliminary  dual  energy  tests  were  also  scanned.  These 
were  placed  side-by-side  with  pieces  of  0.010"  thick  silicon  and  0.015"  thick 
sapphire  taped  to  the  flat  back  surfaces.  A  photograph  is  shown  in  Figure  12. 
This  object  was  intended  to  test  the  ability  of  the  technique^to  distinguish 
materials  that  differ  only  slightly  in  effective  atomic  number.  As  will  be 
seen,  this  phantom,  due  to  an  inadvertent  error  on  the  part  of  the  machinist  who 
made  them,  provided  the  most  dramatic  demonstration  of  Z-sensitivity  of  all  of 
the  objects  studied. 


23 


3UAI  EHEJSY  DIJITOL  RADIOGROPHY 


The  radiographic  scan  data  Is  shown  in  Figure  13.  The  five  steps  on  the 
left  are  one  piece  of  aluminum,  and  the  four  steps  on  the  right  are  a  second 
piece  of  aluminum.  Figure  14  shows  extracted  Compton  and  photoelectric  images. 

A  serious  anomaly  is  immediately  evident.  Since  the  bulk  of  the  material  is 
aluminum,  the  relative  contrast  of  the  two  aluminum  pieces  should  be  similar. 
This  ’s  clearly  not  the  case.  The  aluminum  steps  on  the  right  clearly  have  a 
larger  photoelectric  attenuation  than  the  aluminum  steps  on  the  left.  The 
Compton-to-photoelectric  ratio  for  the  left  stairs  of  aluminum  is  3.7;  for  the 
right  ones,  it  is  2.35.  For  pure  aluminum,  the  expected  ratio  is  4.4.  In  an 
effort  to  interpret  these  results,  the  expected  Compton-photoelectric  ratio  for 
several  common  aluminum  alloys  was  calculated.  On  this  basis,  the  alloy  on  the 
right  was  tentatively  identified  as  either  7075  or  2024.  The  alloy  on  the  left 
was  tentatively  identified  as  6061.  A  spectroscopic  analysis  determined  that 
the  alloys  were  2024  and  6061,  respectively. 

Figure  15  shows  images  with  the  6061  subtracted  and  with  the  2024 
subtracted.  The  left-hand  image  shows  negative  residuals,  i.e.,  regions  with 
higher  apparent  atomic  number  (Z).  Since  the  effective  Z  of  6061  is  nearly  the 
same  as  that  of  silicon  and  higher  than  that  of  sapphire,  we  do  not  expect  to 
see  any  of  the  wafer  fragments  in  this  image.  The  right-hand  image  shows 
positive  residuals  (regions  of  lower  Z  than  Al-2024).  Here  both  the  sapphire 
and  silicon  fragments  can  be  distinguished. 

4.3.3  Air  Force-Supplied  Antenna  Panel.  Figure  16  is  a  photograph  of  a 
corroded  antenna  mount  panel  supplied  by  the  Air  Force  as  a  test  sample.  The 
dual-energy  scans  of  this  sample  are  shown  in  Figure  17.  The  extracted  Compton 
and  photoelectric  images  are  shown  in  Figure  18.  The  panel  was  assumed  to  be 
made  of  the  aluminum  aircraft  alloy  7075.  The  nominal  composition  of  7075  gives 
an  effective  Z  of  15.3.  Figure  19  shows  the  residuals  with  the  Photoelectric 
image  subtracted  from  the  Compton  image  for  effective  Z's  of  15.3  (the  predicted 
value)  and  15.4  (the  value  which  gave  the  best  enhancement  of  the  corrosion). 

The  corrosion  is  visible  as  the  badly-eroded  region  at  the  left  center  of  either 
image.  The  major  portion  of  this  corrosion  is  beneath  the  surface  of  the  panel 
within  what  appears  to  be  stress-induced  delaminations.  The  small  but  definite 
difference  between  the  two  images  gives  a  qualitative  sense  for  both  the 
accuracy  and  the  sensitivity  with  which  the  effective  Z  of  a  material  can  be 
determined. 
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Also  visible  are  some  small  holes  around  the  inner  periphery  of  the  panel 
which  have  been  plugged  with  an  aluminum  alloy  which  has  significantly  lower  Z 
than  the  panel  material. 

4.4  Discussion  of  Results 

The  principle  of  material  identification  and  bulk  material  subtraction  in 
dual  energy  images  has  been  convincingly  demonstrated  in  these  experiments. 
While  it  is,  perhaps,  not  surprising  that  the  Compton-to-photoelectric  ratio 
which  was  used  in  the  generation  of  the  fit  coefficients  should  produce  a  good 
separation  of  the  titanium  from  the  graphite  in  the  images  of  the  calibration 
phantom,  it  is  very  encouraging  that  this  ratio  can  also  be  calculated  quite 
accurately  for  other  materials,  such  as  the  various  aluminum  alloys  studied 
here. 

Several  correctable  deficiencies  are  apparent  in  the  images  obtained. 

First,  improvements  in  both  the  signal-to-noise  ratio  and  the  signal  stability 
would  be  desirable.  Second,  the  thinnest  calibration  phantom  steps  should  be 
thinner  if  relatively  thin  specimens  (such  as  the  antenna  panel)  are  to  be 
optimally  imaged.  Third,  the  x-ray  beam  could  have  profitted  from  additional 
filtration  at  250  keV.  And,  lastly,  a  difference  in  the  x-ray  spot  size  or 
position  manifests  itself  in  the  bright  borders  that  consistently  appear  at 
discontinuities  in  the  image.  Methods  exist  which  can  •'educe  or  eliminate  this 
problem. 
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5.1  SuMMrjf  of  Results 

We  have  successfully  developed  the  techniques  and  generated  the  computer 
algorithms  for  extracting  elemental  information  from  dual-energy  radiographs. 

We  have  demonstrated  both  by  simulation  and  with  real  data  that  the  technique 
can  be  used  to  distinguish  materials  of  different  atomic  number  and  can  enhance 
the  conspicuity  of  minor  components  of  an  object  by  cancelling  the  dominant 
structure  in  the  object.  In  particular,  the  technique  can  be  used  to  greatly 
enhance  the  visibility  of  corroded  material  in  an  image  by  eliminating 
structural  noise  associated  with  more  dominant  material  in  the  component  of 
interest.  Further,  the  technique  is  of  general  validity  and  can  in  principle  be 
adapted  to  any  radiographic  system  capable  of  generating  digital  images  at 
different  effective  energies.  Thus,  it  is  the  general  conclusion  of  this  study, 
that  dual-energy  radiography  can  significantly  enhance  the  performance  of 
real-time  systems  and  deserves  consideration  as  an  adjunct  to  future  and 
existing  radiographic  units. 

5.2  RecoHKndatlons 

Because  the  technique  depends  on  obtaining  accurate  absolute  values  of 
x-ray  transmission,  the  requirements  on  the  data  acquisition  electronics  are 
quite  stringent  but  well  within  the  current  state  of  the  art.  Development  of  a 
linear  detector  array  specially  tailored  for  dual-energy  scans  would  improve  the 
results  obtained  here,  perhaps  significantly.  Thus,  the  development  of  a  linear 
array  tailored  to  and  suitable  for  pursuing  further,  more  quantitative, 
evaluations  of  dual-energy  radiography  is  strongly  urged. 
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APPENDIX 


REVIEW  OF  DUAL-ENERGY  DIGITAL  RADIOGRAPHIC  PRINCIPLES 

There  are  four  basic  types  of  photon  interactions  with  matter:  (1)  inter¬ 
actions  with  atomic  electrons;  (2)  interactions  with  nucleons;  (3)  interactions 
with  the  electric  field  surrounding  nuclei  or  electrons;  and  (4)  interactions 
with  the  meson  field  surrounding  nucleons.  The  effect  of  any  of  these  inter¬ 
actions  can  result  in:  (1)  complete  absorption  of  the  incident  photon;  (2)  elas¬ 
tic  scattering  of  the  incident  photon;  or  (3)  inelastic  scattering  of  the 
incident  photon.  Thus,  in  theory,  there  are  12  different  processes  by  which 
electromagnetic  radiation  can  be  absorbed  or  scattered  by  matter.  In  the  energy 
range  from  a  few  keV  to  many  MeV,  it  turns  out  all  but  a  small  number  of  minor 
effects  are  explainable  in  terms  of  just  three  of  the  twelve  processes.  These 
are  the  photoelectric  effect,  the  Compton  effect,  and  the  pair  production 
effect.  An  excellent  introductory  explanation  of  the  physics  of  these  processes 
can  be  found  in  (7)  and  a  more  advanced  treatment  in  (8). 

An  important  additional  simplification  occurs  when  the  energy  of  interest 
is  well  below  or  well  above  1  MeV.  In  both  instances,  the  physics  of  the 
attenuation  process  simplifies  still  further,  being  dominated  by  just  two  of  the 
three  principle  attenuation  mechanisms.  In  the  former  case,  the  attenuation 
process  can  be  described  to  an  excellent  approximation  as  the  result  of 
photoelectric  and  Compton  interactions;  in  the  latter  case,  as  the  result  of 
Compton  and  pair  production  interactions.  Of  the  two  cases,  only  the  low-energy 
situation  has  been  well  studied,  and  it  is  this  case  which  is  presented  here. 

As  a  result,  we  write  for  present  purposes  the  linear  x-ray  attenuation 
coefficient  as  the  sum  of  two  terms,  one  representing  attenuation  due  to  photo¬ 
electric  interactions,  designated  t,  and  one  representing  attenuation  due  to 
Compton  interactions,  designated  o.  That  is, 

#i(E,S)  =T(E,S)  +a(E,S)  .  (1) 
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The  individual  interactions  in  Eqn.  1  can  be  rewritten  as  follows: 


T(E,S)  »  N{S)  -TaCE.S) 


(2a) 


and 


a(E.S)  =  N(S)  -aaCE.S)  .  (2b) 


where  N(S)  is  the  number  of  atoms  per  unit  volume  at  locality  S  and  and 
are  the  atomic  interaction  cross  sections.  N(S)  equals  the  product  of 
Avogadro's  number,  Nq,  and  the  ratio  of  the  local  material  density,  p(S),  and 
the  local  mass  number,  A(S).  That  is. 


N(S)  » 


Nop(S) 

A(S) 


(3) 


The  atomic  interaction  cross  sections  can  also  be  rewritten  to  explicitly 
exhibit  their  dependencies  on  the  local  atomic  number,  Z(S),  and  the  energy  of 
the  incident  photon,  E.  From  references  (7,8),  the  atomic  cross  sections  can  be 
represented  with  sufficient  accuracy  for  present  purposes  as: 


Ta(E,S)  =  K^z4(s)  f^(E) 


(4a) 


and 


Oa(E,S)  =  Kj,Z(S)  f(,(E)  ,  (4b) 

where  fj  and  fg  carry  the  energy  dependencies  of  the  interactions  and  and 
are  dimension-carrying  constants  independent  of  photon  energy  or  material  compo¬ 
sition.  Thus, 

t(E,S)  »  •  S-Z'^(S)  fT(E)  (5a) 

A(S) 


and 


HgPiS) 


.  KjyZ(S)  f^(E) 


a(E,S)  = 


(5b) 


The  crucial  feature  of  Eqn.  5  Is  the  fact  that  the  spatial  and  energy 
dependencies  (S  and  E)  separate.  That  is,  Eqn.  5  can  be  rewritten  as 


t(E,S)  =  a^(S)  f^(E) 

(6a) 

a(E,S)  =  a^(S)  f^(E) 

(6b) 

where  a^  and  a^  carry  strictly  the  spatial  dependencies  shown  in  Eqn.  6  and  fj 
and  ffj  carry  strictly  the  energy  dependencies.  With  the  aid  of  Eqn.  5,  Eqn.  6 
can  also  be  expressed  in  the  following  manner: 

T(E.S)ocp(S)Z3(S)  •  f^(E)  (6c) 

and 

a(E,S)«p(S)  .  yE)  .  (6d) 

where  Eqn.  6c  and  6d  makes  note  of  the  fact  that  the  ratio  Z/A  is  approximately 
constant  over  the  entire  periodic  table  (except  for  hydrogen).  Consequently, 
the  photoelectric  attenuation  is  roughly  proportional  to  pZ^  and  the  Compton 
attenuation  to  p. 

The  essential  result  of  this  discussion  can  now  be  derived.  In  radiog¬ 
raphy,  the  observed  x-ray  intensity,  1(E),  can  be  related  to  the  incident 
intensity,  Io(E),  by  the  well-known  law  of  radiation  physics: 

1(E)  =  Io(E)  exp  .S)dS]  ,  (7) 

where  Jp/u (E,S)dS  represents  the  line  integral  of  the  linear  attenuation  co¬ 
efficient  of  the  object,  fi(E,S),  along  the  ray  path  r  defined  by  the  particular 
projection  of  the  object  being  taken.  For  a  monochromatic  flux,  the  line 
integral  of  along  the  given  ray  path  can  be  readily  extracted  by  forming 

the  negative  natural  logarithm  of  the  ratio  of  I(E)/Io(E), 

L(E)  =  -  in(I/Io)  =/pM(E,S)dS  .  (8) 

When  the  incident  spectrum  is  polychromatic,  the  prescription  given  by  Eqn.  8 
for  extracting  the  value  of  /judS  is  still  valid  provided  E  is  interpreted  as 
the  effective  energy  of  the  x-ray  bean.  Since  the  effective  energy  varies  as 
the  ray  path  through  the  object  changes,  the  calibration  procedure  becomes  more 
complicated,  the  overall  validity  of  the  method  is  not  affected. 


^3  "Jl  'J»  "Ji  ' 
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utilizing  the  results  of  Eqns.  1  and  6,  Eqn.  8  can  be  rewritten  in  the 
following  more  convenient  form: 


*/j,^(E,S)dS  =/r(T(E,S)  +a(E.S))dS 

(9a) 

=  f^(E)/j,a^(S)dS  +  f^(E)/ra^(S)dS 

(9b) 

*  ciL^  +  , 

(9c) 

where  in  Eqn.  9,  f^CE)  and  f0(E)  have  been  replaced  by  constants,  ci  and  C2,  to 
emphasize  that  they  depend  solely  on  the  effective  energy  of  the  incident 
radiation  and  L  and  L  are  line  integrals  along  the  same  ray  path  and  depend 
solely  on  the  material  properties  of  the  object. 

The  implication  is  as  startling  as  it  is  simple;  if  scans  are  performed  at 
two  different  effective  energies,  E  and  E',  then  l  and  \?  can  be  extracted 
separately,  point-by-point,  by  solving  the  following  pair  of  simultaneous 
equations: 

L(E)  =  Cl  l’’  +  C2  (10a) 

L(E')  =  ci'C  +  C2'L‘^  .  (lOb) 

If  the  data  sets  from  two  identical  x-ray  projections'*",  iL(E)[  and 
|l(E')|,  were  noiseless,  then  any  two  different  effective  energies  would 
suffice;  but  of  course,  in  any  non-mathematical  context  they  cannot  be.  This 
implies  that  some  care  must  be  taken  to  engineer  (to  the  extent  possible)  a 
situation  that  will  produce  the  greatest  degree  of  linear  independence  between 
the  above  equations.  Stated  another  way,  the  solution  to  Eqn.  10  will  depend  on 
the  statistical  accuracy  of  the  measurements,  L(E)  and  L(E'),  and  the  nature  of 
the  determinant  of  the  coefficients  on  the  right  hand  side  of  the  equation.  It 
is  assumed  that  the  measurements  can  be  engineered  so  that  ciC2'-ci'c2  is  well 
conditioned  and  non-zero,  and  that  exposures  are  sufficient  to  obtain 
well-behaved  data  sets.  In  practice,  there  is  sufficient  flexibility  that 
these  are  not  severe  constraints.  Thus,  in  the  following  development,  it  will  be 
supposed  that  the  two  effective  energies  are  well  separated  and  that  the 
measurements  are  statistically  accurate,  two  demands  which  are,  at  least 
intuitively  if  not  quantitatively,  well  defined. 


^The  notation  adopted  here  uses  brackets  to  indicate  a  collection  of  ordered 
points  which  constitute  an  image. 
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Under  these  conditions,  Eqn.  10  can  be  solved  straight-forwardly, 
point-by-point,  over  the  whole  image  with  the  result: 

|L(E)},  {L(E')[=:^  {J’},  {L*"}  .  (11) 

In  words,  Eqn.  11  states  that  the  two  new  images,  {L^|  and  ,  can  be 

generated  from  the  dual-energy  images,  {L(E)[  and  |L(E')|,  to  obtain  radio¬ 
graphs  proportional  to  the  material-dependent  portion’s  of  the  projected  object. 
It  also  states  further  that  since  a^op  (S)z3{S)  and  (S),  the  so-called 

dual-energy  images  provide  a  powerful  method  of  sampling  the  density  and  atomic 
number  within  an  object  of  interest.  It  is  this  potential  access  to  two 
fundamental  material  attributes  that  make  dual-energy  radiography  such  an 
attractive  technique  for  corrosion  detection. 

To  see  how  this  would  work  in  practice,  consider  the  following  important 
extension  of  the  above  derivation.  Eqn.  11  states  that  a  high-  and  a  low-energy 
images,  |L|  and  jL'|  ,  can  be  decoupled  to  produce  two,  more  fundamental 
images,  one  proportional  to  the  density  of  the  projected  material  and  one 
proportional  to  the  product  of  the  density  and  the  cube  of  the  atomic  number  of 
the  projected  material.  However,  and  this  is  the  key  point,  these  two  derived 
images,  |l  |  and  ,  can  in  turn  be  manipulated  to  produce  an  almost  end¬ 

less  series  of  new  images  which,  though  less  physical  in  nature,  are 
nevertheless  considerably  more  practical  from  the  standpoint  of  corrosion 
detection. 

The  basic  procedure  is  to  form  linear  combinations  of  {L^j  and  {L*^!  in 
such  a  way  as  to  achieve  some  useful  advantage.  Since  this  is  done  point-by¬ 
point,  the  brackets  can  be  dropped  for  convenience  without  loss  of  generality  in 
the  following  discussion.  First,  define  a  general  linear  combination  of  the 
photoelectric  and  the  Compton  image  as  follows: 

R  =  kL^  +  k’L*'  ,  (12) 

where  k  and  k'  are  constants  and  R  stands  for  a  point  in  the  resultant  (or 
residual)  image.  It  can  be  seen  immediately  by  comparing  Eqn.  12  with  Eqn.  9c 
that  if  k  is  taken  equal  to  cj  and  k'  to  C2,  R  just  becomes  one  of  the  original 
images.  It  may  appear  as  though  we  have  just  reversed  our  earlier  efforts;  but 
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properly  viewed,  this  result  is  far  from  trivial.  It  implies  that  by  taking 
physically  meaningful  values  for  k  and  k‘,  we  can  produce  radiographs  at  any 
effective  energy  we  wish  from  our  two  dual-energy  images.  This  might  be  very 
important  in  scanning  for  low-contrast  detail  in  some  applications. 


On  the  other  hand,  other  more  general  linear  combination  can  also  be 
formed.  Perhaps  the  most  powerful  of  these  approaches  for  corrosion  detection 
is  the  one  that  seeks  to  cancel  a  designated  material.  Imagine  for  concreteness 
that  an  object  consists  predominantly  of  aluminum.  In  this  case,  the  line 
integrals  shown  in  Eqn.  9  can  be  written  as  a  term  integrated  over  the  projected 
aluminum  material  and  a  term  integrated  over  whatever  residual  material  is  also 
present.  Carrying  the  separation  forward,  Eqn.  12  can  be  rewritten  as: 

R  =  k(L‘^(Al)  +  L^(res))  +  k'(L‘^(Al)  +  L^^Cres))  (13a) 

=  [kL^ (Al)  +  k'L^^CAl)]  +  [kL^  (res)  +  k'L^^Cres)]  ,  (13b) 


where  L(A1)  and  L(res)  stand  for  the  appropriate  integrals  over  the  projected 
amounts  of  aluminum  and  residual  materials,  respectively.  By  examination,  it 
can  be  seen  that  by  taking: 


ill.  =  .  L^(Al) 
k  L‘^(A1) 


(14) 


the  aluminum  contribution  to  any  point  R  can  be  made  to  vanish.  Thus,  by 
performing  the  subtraction  on  a  point-by-point  basis  over  the  whole  image,  the 
residual  (non- aluminum)  material  in  the  object,  R,  can  be  viewed  with  all 
aluminum  detail  "magically"  cancelled. 

This  is  obviously  a  very  significant  and  powerful  result.  For  one  thing, 
it  implies  that  dual-energy  techniques  can  be  used  to  remove  dominant  structure 
in  an  image  and  reveal  low-contrast  detail  that  might  not  otherwise  be  observed 
due  to  structural  noise.  For  another,  it  suggests  that  this  method  can  also  be 
used  to  help  identify  foreign  materials  within  an  object.  To  see  this  more 
clearly,  consider  Eqn.  14  again.  Since  is  proportional  to  pZ^  and  is 
proportional  to  p,  it  follows  that 

JilocZ^  (15) 

k 
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In  general,  one  would  not  expect  to  have  to  identify  the  internal  detail  in 
an  object.  That  would  be  known  a  priori  from  assembly  drawings.  However,  if  a 
foreign  substance  or  unexpected  object  were  to  appear  in  a  radiograph,  the 
operator  could  adjust  k  and  k*  (perhaps  automatically  with  a  joystick  or  track¬ 
ball)  until  the  item  of  interest  exactly  disappears  from  the  residual  image. 

The  atomic  number  of  the  item  could  then  be  inferred  from  the  ratio  of  k  and  k' 
using  Eqn.  15.  It  is  conceivable  that  this  procedure  could  prove  useful  for 
assessing  corrosion  progression  since,  in  some  circumstances,  the  severity  of 
the  corrosion  present  might  be  proportional  to  the  relative  amounts  of  oxide, 
salt,  and/or  hydride  formed.  That  is,  dynamics  of  the  situation  might  be  such 
that  shifts  in  the  average  atomic  number  of  the  corroded  material  could  be 
correlated  with  the  degree  of  progression  of  the  damage. 
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